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■ We examine the gravitational collapse of an infinite cylindrical distribution of time like dust. In 

5^ I order to simplify the calculation we make an assumption that the axial and azimuthal metric 

functions are equal. It is shown that the resulting solution describes homogeneous collapse. We 
show that the interior metric can be matched to a time dependent exterior. We also discuss the 
^ ■ nature of the singularity in the matter region and show that it is covered. 

a^ 

00 

CN ■ 1 Introduction 
(N 

' While there have been many analytical studies of spherical gravitational collapse, there have been 

only few investigations of non-spherical collapse, the obvious reason being the complexity of Ein- 
^ ■ stein equations. Thus most of our present understanding of cosmic censorship and the nature of 

bX)' singularities stems from studies of spherical collapse. Given the impact that the results of these 

studies have had on our understanding of the end states of classical collapse, it is of interest to ask 
^ , what may emerge from studies of gravitational collapse with different topologies. 

a I A very useful study of the collapse of an infinite dust cylinder was initiated by Thorne . This 

is described by the metric 



ds' = e(^-^) {dt^ - dr^) - e^dz^ - a'e-Pd(t>\ (1) 

The metric functions depend on t and r. This form of the metric can be inferred by first writing 
down the static Weyl axisymmetric metric for an infinite line Newtonian source 

in which the metric functions depend only on r and z. The transformation t ^ iz, z ^ it in this 
metric then yields the metric It is invariant under Lorentz transformations in the {t,r) plane 
and can be used to represent the collapse of a cylindrical dust shell. Apostolatos and Thorne P] 
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used the above metric to show that rotation can halt the collapse of the cylinder. The formation 
of a naked singularity for the infinite cylindrical shell collapse was established by Echiverria [H]. 
Naked singularities for counter-rotating dust shells were shown to occur, by Goncalves and Jhingan 
[3] and by Nolan 5 , and an exact dynamical solution was given by Pereira and Wang The 
collapse of infinite cylindrical dust clouds has been studied also by Chiba [7j, by Senovilla and Vera 
[H] and by Bondi 0. 

Although the collapse of an infinite cylinder may not constitute a "realistic" collapse scenario, 
it can simulate the collapse of a finite "bar" or spindle-like matter distribution very near the central 
regions of the spindle. In this sense the collapse of an infinite cylinder is of astrophysical interest. 
More importantly, however, the problem probes the structure of the general theory of relativity. 
Our analysis in the present paper is of an infinite cylindrical cloud of time-like dust and is motivated 
by what is known for spherical time- like dust collapse, i.e., the Tolman-Bondi spacetime. We make 
one assumption, namely that the azimuthal and axial metric components are equal to each other. 
As we show, this assumption implies that the solution is homogeneous. Also, the resulting system 
is a 2-d system in which the metric depends on a radial coordinate, r and on time, t. 



2 Collapse of an infinite dust cylinder 

Consider an infinite cylindrical cloud of pressure-less, inhomogeneous, time-like dust described by 
the stress tensor 

T^u = e{x)u^Uy, (3) 

where = —1. We first set up the metric for the interior of this infinite dust cylinder, using 
comoving coordinates (t, r, z, (p). Assuming that z represents the axis of symmetry of the cylinder, 
the metric functions depend only on t and r. The dust cloud is assumed to extend up to some radius 
r^; matching to a vacuum exterior will be discussed in Section 5. As a consequence of axisymmetry 
the functions gtrj,, gr<j) and gz<j, all vanish. Moreover, for an infinite cylinder, invariance under z —z 
implies that gtz and grz also vanish. We can use the freedom of coordinate transformations from 
the pair (t, r) to a new pair (t', r') to set two functions to zero: these are gtr and the radial velocity 
u*". As a result of the latter condition, r is determined to be a comoving coordinate. Since the 
matter is dust, the metric can be comoving and synchronous, so that (700 = — 1- The metric for an 
infinite dust cylinder can hence be written in terms of three unknown functions of two variables as 

ds^ = -df + L^{t,r)dr'^ + M'^{t,r)dz^ + B'^{t,r)d(j)'^ . (4) 

The energy-density e{t, r) being the only non-zero component of the energy-momentum tensor, 
Einstein's equations for this spacetime are (dot and prime denote derivatives w.r.t. t and r, respec- 
tively) : 

_ M B_ MB _ B'M' 

B^ MB~ BML^' (5j 

B L BL B"L-B'L' 
= B^L^B-L= BL^ ' 
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M L ML M"L-M'L' 
^^^-^'■M^L^ML- ML3 ' 

M' B' L (W B'\ 

M B L\M B ) ' ^ ' 

_ B"L - B'V M"L - M'V B'M' BL ML B M _ 

Gtt-e{t,r): __ + __ + __ + ___£. (9) 

The conservation of the energy-momentum tensor gives the relation 

where ip{r) is an integration function, to be determined by the initial data. We also note that 
adding equations © , © and ((Tj) , and using @ and ifTU]) gives the relation 

(u) 

M B L LMB ^ ' 

Instead of working with the five Einstein equations ©-(dl) we will work with the equivalent set 
given by the equations ®, ©, (CHI), CU) and the difference (©-©• 

We now assume that B{t,r) = rM{t,r). For a finite axisymmetric object like a spheroid this 
would probably not be allowed by Einstein equations but, as we see below, the equations are self- 
consistent when this assumption is made for an infinite cylinder. The study of this system may be 
thought of as a prelude to examining the most general case, when B and M are not related. The 
physical meaning of the assumption is that the object shrinks at the same rate along the radial 
direction and the axis, so that its 'prolateness' or 'oblateness' does not change with time. More 
importantly, we show that this assumption implies that the matter distribution is homogeneous. It 
is not the most general homogeneous solution, but a special case. 

With the above assumption the metric (|1J) becomes 

ds"^ = -df + L^it, r)dr^ + M^{t, r) [dz^ + r'^dcj?] (12) 

and Einstein's equations for this metric are as follows: eqn. @ becomes 

2B B^ 1 B' 



while eqn. (|TT1) reduces to 



B ^ B^ L'^Kb'^ Br' ' ^^^^ 

2B L Mr) , , 

B ^ L rLM'^ ^ ' 

These two are the dynamical equations, while the remaining three equations are constraints. Eqn. 
(jHI now becomes 

2M' M L f2M' r 



2M' L' 



M rM L\M r . 
whereas the difference © — gives 

M' T! 

(16) 



M L 
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and the conservation equation H1U|) now is 

Eqns. H13p -()17 p are the Einstein equations for the metric (|12() . Eqn. 1)16(1 can be solved to give 

L{t,r) = h(t)M^{t,r) (18) 

where h{t) is a arbitrary function of time. Eqn. ()15|) imphes that 

(V^M)' = 5(r)L(t,r) (19) 

where (7(r) is an integration function. Using (jl9j) and ()18() in the first dynamical equation (^Sj) 
gives 

2iji3 + ij2 = r52(r)-^^ (20) 
while the second dynamical equation becomes 

i?/. + 4i?/.-— ^ + ^--^. (21) 

3 Exact solution for a class of homogeneous dust collapse 

We choose the scaling 5(0, r) = r at time t = 0, when collapse is assumed to begin. This implies 
M(0, r) = 1. We also choose /i(0) = 1, which from (|18() implies L(0,r) = 1. From Eqn. p9() . 
written at t = 0, we get the relation g{r) = 1/2-y/r. Then, solving (fT9|l yields 

l^ = h{t)+K{t)V^. (22) 

In this equation, h(t),K{t) are functions of time only. Now the differential equations ()20() and 
(PT|) should be valid at all r (including r = 0). Substituting B = rM in equation gives 

2M + - ^ 



4r2 Ar'^h'^M'^ 
h^-{h + K^f 



(23) 



4^2 /j2 

The l/r^ term in the denominator makes the equation invalid for r = 0. This can be avoided by 
setting K{t) = 0, which implies h = 1/M and hence that M is a function only of time. Similarly 
equation (|21|1 gives 

M/l + 4M/t = -^+ , I , (24) 

The sum of the last two terms vanishes when we set K{t) = 0. There is a 1/r in the denominator 
of the coefficient of ip{r). Since the left hand side is independent of r, we are constrained to choose 
V'(r) = kir, where ki is a constant. Substituting for ip in equation (|17() gives 

e(i,0 = /^ = 2/ci/i3. (25) 
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Since h is a function of time only, this means that the density distribution is homogeneous. Exphcit 
solution for the metric can be obtained by solving the equation 



2MM + = (26) 

which gives 

M = {cit + C2)2/^ (27) 

Substituting the value of M in equation pi]) gives the relation between ki and ci; ki = 2cf/3. 
Also, with h{0) = 1, and assuming ci to be negative, we can write 

M = c/\to-tf/^ (28) 

where Iq = l/|ci|. The metric can now be written as 

ds'^ = -df + k^ito - t)^/^{dr^ + dz^ + r^d(t>^) (29) 

where A; is a constant. At time t = to there occurs a curvature singularity. All the shells become 
singular at the same time Iq. This is similar to the spherical Oppenheimer-Snyder dust collapse. 



4 Comparison with spherical dust collapse 



Let us compare this system of equations with those for spherical dust collapse, described by the 
Tolman-Bondi spacetime: 

ds'^ = -dt^ + {t, r)dr^ + {t, r) [d9^ + sin^ Odcp'^] . (30) 

The Einstein equations for this metric are: 

2ij i?2 1 /i?'2 X 
G..=0:- + ^ = ^(^-lj, (31) 

L R L R R" L — R' L' 
Gee = G^^ = 0:- + - + --= , (32) 

Gtr = 0:^ = ^, (33) 

R"L-R'L' R^ LR 1 _ 

^" - ' ^ RL^ 2m^ ^W^^LR^W^-' ^^^^ 



The conservation equation is 



e(i,0 = ^- (35) 



We also note that adding equations (|HT|) and and using (jHH) gives 



2R L Mr) , , 

If + 1 = -IT? ■ P"' 
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These equations should be compared with those for the cyUnder: Eqn. (|5T|) should be compared 
with Eqn. (fT^ . Eqn. I\32\i should be compared with Eqn. ©, Eqn. should be compared with 
Eqn. (|15|) . Eqn. (|34|) should be compared with Eqn. Q, Eqn. I|35|) should be compared with Eqn. 
(|17l) . and Eqn. (|36|) should be compared with Eqn. (|14j) 

Thus, although there is similarity in the forms of the metrics ()1'2() and (|3U() . the cylindrical 
metric ()12[) describes homogeneous collapse, whereas the spherical metric ()3Up describes general 
inhomogeneous collapse. This is because of the additional Einstein equation (|16|) which is not there 
in the spherical case. 



5 Matching the interior cyUndrical metric to a static exterior 

The metric H'29() is valid in the interior of a collapsing dust cloud. The exterior is vacuum. The 
above metric will be a solution of Einstein equations if it can be successfully matched to a static 
or a time dependent exterior. It can be easily guessed that the interior metric might not match 
to a static exterior, since the quadrupole tensor for the matter distribution is non-zero, and the 
collapsing system emits gravitational waves. 

To verify explicitly the impossibility of static exterior vacuum solution, we assume a static 
exterior metric of the form 

ds"^ = -e^dT^ + e^dR^ + e'^'dZ^ + R^d<S>^ (37) 

where A,B,C are functions of R. Only the diagonal components of the Einstein tensor survive. 

e^~^{-2C' -RC''^ + B'{2 + RC')-2RC") , , 

Goo = ^ (38) 

2a^A!{2^Ra) 
Gu = ^ (39) 

e^-^{RA''^ -2B' + A'{2- RB') + 2RA") 
G22 ^ (40) 

e-^R^{A'^ - B'C + C'^ + A'{-B' + C') + 2A" + 2C") . 
G33 = ^ (41) 

All the components equal zero in vacuum. Solving the Einstein equations gives 

G = /iln(^-|-) +/2. (42) 



-2/1 In 



'1 1" # 

B^^.k. (44) 
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Rq is a length scale which, as seen below, gets fixed by the matching. We can get rid of I2 and 
^3 by rescaling T and Z. So there are two unknown constants, (see also [S]). 

In the interior let the outermost shell be labelled by r = r<j. The coefficient of is interpreted 
as the square of the radius. The expression rsk{tQ — t)"^/^ represents the evolution of the radius of 
the outer boundary of the collapsing dust with time. Further, it can be assumed that $ and Z 
represent the same coordinates in both interior and exterior, i.e. ^ = (f) and Z = z. Consider the 
hypersurface r = rg. On this hypersurface, the interior metric is given by 

ds^ = -dt^ + k\to - tY'^idr"^ + dz^ + rldcf^). (45) 

The exterior is given by (jHTj). 

The metric coefficients of (j) and z should be same in both the interior and exterior. Equating 
the coefficient of d</>^ we get 

Rl = rlk\to-tf^ (46) 

The subscript s in Rg is to indicate the hypersurface in the exterior coordinates. Similarly equating 
the coefficients of dZ"^ for both interior and exterior gives 

which implies li = 2 and Rq = r^. So the exterior metric on the hypersurface is 

ds^ = -^dT^ + ^dR^ + f ^ y + r^d^^ (48) 
Rs Ts \rs J 

where c is a constant to be determined by matching the second fundamental form. Matching the 
remaining components of the metric on the hypersurface yields 

-dt' = -I^dT^ + 2^dR^, (49) 

Rs rs 

This can be used to obtain the relation between the interior and exterior time. The second funda- 
mental form JU] is given by 

$ = {-n^.ydx^'dx'') (50) 

and the extrinsic curvature is given by 

K^u = n^.y (51) 

where is the unit normal to the hypersurface. For the interior metric given by ()45() the hyper- 
surface is given by r = r<j. Calculating K'^ yields 

where Rs is given by ()46() . Similarly for the exterior metric K"^ can be computed. The normal to 
the hypersurface Rg can be obtained by differentiating (|46|) . 

(2rsVk(to-t)-^/^^)dT , , 

dR + ^ ^ J ^ — = 0. (53) 
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Using the above equation and evaluating K'^ gives 

< = t (54) 

where is the contravariant component in the R direction. By equating K'jj^ for both interior and 
exterior we obtain the value for c, 

^ (55) 



Ts 9 

c was initially assumed to be a constant. In equation ()55() c is dependent on Rg which is given by 
H46|) . So c is not a constant as assumed but is time dependent. Hence there is no static exterior for 
the metric given by (|45|1 . 



6 Matching with time dependent exterior 

In this section it will be shown that the interior metric ()45() can be matched to a time dependent 
exterior. Let us assume an exterior of the form, 

ds^ = e(«-^)(dT2 - dR^) - e^'dz^ - a^e'^^^, (56) 

where G, P, a are functions of R, T. The Einstein equations are 

a'G-aP'P + G'a-2a' ^ 

Gtr = = U (5»j 

Za 

O,«=-"'P" + ^')+^'Q'°'-^° + Q°)=0 (59) 

4a 



e-^+2-P(a(P'2 + 2G"-4P"-p2_2G' + 4P)+4(-PV + a" + Pd-a)) ^ 
Gzz = = (60) 

e-ga^(P-^ + 2G---p2-2G) 

tr00 = ^ = U (61) 

Here prime and dot denote partial derivative with respect to R and T respectively. In vacuum all 
the above expressions equal zero. Subtracting equations H57() and (|59() gives 

q" - Q = (62) 

Equation (|HT|) gives 

p2 _ p/2 

G" - G = (63) 
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Substituting equations and H63() in equation ((Sn|l gives 

P'a' - Pa 

P"-P+ = (64) 

a 

The remaining two equations can be written as 

Ga - P'Pa + G'a - 2a = (65) 

-a(P'2 + P2) + 2{G'a' - 2a" + Go) = 0. (66) 

The equations (|62|). (|63j) . I|64|) are second order partial differential equations in T and R. The 
outer boundary of the collapsing matter (hypersurface r = r^) is a priori unknown in the exterior 
coordinates. It is of the form /i(r, R) = const. To solve completely the partial differential equation 
(|62|) (second order in both T and i?), we need to know the value on the hypersurface and also the 
first derivative in the direction perpendicular to the hypersurface. We assume that Z and $ 
represent the same coordinates in both the interior and exterior. So on the hypersurface the metric 
components have to match. Equating coefficient of for both metrics we get 

^ = rlk\t,-tfl^ (67) 
Similarly equating coefficient of dZ"^ we get. 

e^ = fe2(to-t)'/' (68) 

Multiplying both we get 



a 



k^sito-t)"/^ (69) 



To solve the equations of the type a" — a = 0, we already know the value on the hypersurface 
given by ()69() . We need one more arbitrary function to provide information about the value of the 
derivative in the perpendicular direction. Similarly we know the value of P on the hypersurface 
(|68|) . To solve the equation (|64j) we need one more arbitrary function. For solving equation (jHSJ we 
need two arbitrary functions since we do not know the value of G on the hypersurface. These four 
arbitrary functions along with the initial hypersurface fi{T, R) makes five arbitrary functions. The 
matching will be possible if the five arbitrary functions are enough to satisfy all the constraints 
which are required during matching. The relation between the interior time and the exterior 
coordinates can be obtained by the equation similar to equation H49() 

Equations (|65|) and H66j) act as constraints and can be used to eliminate two of the five arbitrary 
functions. Now matching the second fundamental form jJOl 

$ = Kudt^ + K,,dz^ + K^^dcl)^ (71) 

gives at most three independent constraints. There are three arbitrary functions left which can 
take care of the constraints offered by the second fundamental form. So in principle the interior 
metric ()45|) can be matched to the time dependent exterior given by H56() . 
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7 Nature of the singularity 



The question of interest is to find out whether the curvature singularity resulting in the matter 
interior from the gravitational collapse of the dust cloud is naked or covered. This question can be 
addressed by looking at the null geodesic expansion Q for the radial null geodesies. The non-zero 
components of the tangent to the null geodesic are K^{t,r) and K^{t,r). The calculation of for 
the cylindrical metric 1)12(1 proceeds exactly as for the spherical case discussed in and we get 
the result 



2K^- 
Z 



1 + 



B 



Vrg{r) 



2K'- 
Z L 



1 + 2/i(0)^ 



r 



1 + 2rM 



(72) 



where Z = Bj ^ = y^M(t). Z' jZ is positive, because of equation (|T9|) . and K'^ is positive because 
we are examining outgoing geodesies. 



and can be evaluated explicitly from the geodesic equations 



~dk 



M 
M 



dK* 
~dk 



-2M {K 



r\2 



(73) 



Here, k is the affine parameter and we have used the relation dt/dr = M{t) = /K"^ for outgoing 
geodesies, which follows from using ds^ = for radial null geodesies. 

The explicit solution is however not necessary for our purpose. As the singularity is approached, 
i.e as t tends to t^, M tends to -co, and hence so does the geodesic expansion 9, showing that the 
singularity is covered. The expression for 9 is in fact identical to that for the spherical marginally 
bound homogeneous dust solution [Oppenheimer-Snyder] for which the metric is 



ds^ 



-dt^ + S'^{t, r) [dr"^ + r'^d9'^ + sin^ 



(74) 



and the scale factor S{t) goes as (to — t)^/^. Comparison of this metric with the cylindrical ho- 
mogeneous metric (|29() makes it clear that the geodesic expansion 9 behaves identically in the two 
cases. 

It has earlier been demonstrated that the collapse of an infinite dust cylinder results in a 
naked singularity, in accord with the hoop conjecture. How is our result above consistent with these 
earlier findings? A possible explanation is that there is a singular Cauchy horizon in the exterior, 
as discussed recently by Tod and Mena jT^. It should be noted that our analysis here does not 
provide information about the formation of a singularity, or otherwise, in the exterior spacetime. 



8 Conclusion 

In this paper we have considered the collapse of an infinite cylindrical cloud of time-like homoge- 
neous and pressureless dust subject to the condition that the axial and azimuthal metric functions 
are equal. We have explicitly shown that the resulting solution is homogeneous and that the col- 
lapse terminates in a curvature singularity. We have shown that there exists a vacuum exterior and 
it can be demonstrated that the singularity forming in the matter region is covered. 
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